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Abstract 

Amphiphiles are molecules which have both hydrophilic and hydropho- 
bic parts. In water- and/or oil-like solvent, they self-assemble into ex- 
tended sheet-like structures due to the hydrophobic effect. The free energy 
of an amphiphilic system can be written as a functional of its interfacial 
geometry, and phase diagrams can be calculated by comparing the free en- 
ergies following from different geometries. Here we focus on bicontinuous 
structures, where one highly convoluted interface spans the whole sample 
and thereby divides it into two separate labyrinths. The main models for 
surfaces of this class are triply periodic minimal surfaces, their constant 
mean curvature and parallel surface companions, and random surfaces. 
We discuss the geometrical properties of each of these types of surfaces 
and how they translate into the experimentally observed phase behavior 
of amphiphilic systems. 

1 Surfaces in self-assembling amphiphilic sys- 
tems 

The subject of this article is the theoretical description of amphiphilic systems, 
which are one example of soft condensed matter. Soft matter is material that 
has a typical energy scale of fcsT « 4 x 10^^^ J, where ks is the Boltzmann 
constant and T « SOO-ftT is room temperature. For such material, perturbations 
arising from the thermally activated movement of its molecular components are 
sufficient to induce configurational changes, and entropy is at least equally im- 
portant as energy in determining its material properties. The understanding 
of the underlying mechanisms is essential for the application of many technolo- 
gies in everyday life, including colloidal dispersions (paints, inks, food, creams. 



lotions), foams (beverages), liquid crystals (displays), polyelectrolyte gels (dia- 
pers) and soaps (washing and cleaning). Moreover, the term soft matter also 
includes most biomaterials, for example blood or cartilage. One particularly 
important biomatcrial is the biom,emhrane, that is the protein carrying lipid 
bilayer which surrounds each cell and its organelles. A good model system for 
the structural properties of biomembranes are lipid bilayers, which form sponta- 
neously in mixtures of water and lipids and which are one of the main subjects 
of this paper. 

Soft matter systems are very often characterized by competing interactions 
(including cntropic ones) , which lead to structure formation on the length scale 
between tens and hundreds of nanometers (1 nm = 10~^ m). Since this length 
scale is not accessible by optical lithography, self-assembly in soft matter sys- 
tems is one of the main concepts of nanoscience. With structure formation 
being so prominent in soft matter systems, their theoretical description often 
centers around their spatial structure in three dimensions. For a rough clas- 
sification of the different approaches used, it is convenient to use an analogy 
to the field of random geometries, and to distinguish between point, fiber and 
surface processes. Examples of systems which are suited for approaches along 
these lines are colloidal dispersions, polymer networks and sheet-like structures 
in self- assembling amphiphilic systems, respectively. Here we will treat the lat- 
ter case, but we will also show that surface dominated system can be related 
to Gibbs distributions for scalar fields, namely through the use of isosurface 
constructions. 

Amphiphilic systems are solutions of amphiphiles in suitable solvent. Am- 
phiphiles are molecules which consist of two parts, one being hydrophilic (water- 
like, called the head) and one being hydrophobic (oil-like, called the tail). Well 
known classes of amphiphiles are tensides (used for washing and cleaning pur- 
poses) and lipids (the basic components of biomembranes) . Fig. 1 shows several 
examples for which phase diagrams are discussed below. Tensides are often 
called surfactants, due to their surf ace activity at interfaces between water-like 
and oil- like phases. Hydrophilic and hydrophobic molecules demix due to the 
hydrophobic effect; oil-like molecules are expelled from a region of water-like 
molecules since they disturb their network of hydrogen bonds. Hydrophilic and 
hydrophobic parts of an amphiphilic molecule cannot demix due to the covalcnt 
linkage between them, but the amphiphiles can self-assemble in such a way as 
to shield its hydrophilic parts from its hydrophobic ones and vice versa. Am- 
phiphilic systems have many similarities with diblock copolymer systems, for 
which the term microphase separation is used to denote the molecular tendency 
for structure formation (see the contribution by Robert Magerle). However, 
in contrast to diblock copolymers, amphiphiles are small molecules, thus the 
entropy of their molecular configurations plays only a minor role in determin- 
ing the overall structure. Moreover, whereas in diblock copolymer systems one 
component is sufficient to obtain stable mesophases, in amphiphilic systems usu- 
ally the presence of an aqueous solvent is necessary to obtain well-pronounced 
structure formation. Both diblock copolymer and amphiphilic systems show 
structure formation on the nanometer scale, but the length scale is set by dif- 
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Figure 1: Schematic representation of different amphiphilic molecules. Hy- 
drophobic tails are to the left, hydrophilic heads to the right, (a) Pentaethylene 
glycol dodecyl ether, Ci2H2^{OC H2C H2)hOH or in short C12E5, is a small sur- 
factant, (b) Laurie acid, CH^{CH2)wCOOH or in short LA. Fatty acids are 
the simplest amphiphiles and can be dissolved to high amounts in phospholipid 
bilayers. (c) Dilauroyl phosphatidylcholine, or in short DLPC. Phosphatidyl- 
cholines have two hydrocarbon tails and a zwitterionic headgroup. They are 
abundant in animal cells and the best studied model system for biological lipids, 
(d) Monoolein, a monoacyl glycerol with one (unsaturated) hydrocarbon tail. 



3 



Large figure missing ... 



Figure 2: Self-assembly of amphiphiles: (a) in binary systems, amphiphiles self- 
assemble into bilayers, micelles and vesicles in order to shield the hydrophobic 
tails from the aqueous solvent, (b) In ternary systems, the oil- like solvent swells 
the hydrophobic regions. If the oil-like solvent is the majority component, in- 
verse structures occur. In ternary systems, all interfaces are amphiphilic mono- 
layers. 



ferent control parameters: for diblock copolymers and amphiphiles, these are 
polymer size and solvent concentrations, respectively. Since the solvent in an 
amphiphilic system usually has no special physical properties by itself, its prop- 
erties are mostly determined by its interfaces. This stands in marked contrast 
to the case of diblock copolymer systems, where the entropy of chain configura- 
tions in the regions away from the interfaces cannot be neglected in a physical 
description (for a review see [71] and references therein). As we will discuss in 
more detail below, interface descriptions have been very successful in describing 
the properties of amphiphilic systems (for reviews see [63, 40, 64, 104, 92]). 

Amphiphilic systems can be classified according to the solvent used. A mix- 
ture of water and amphiphile is called a binary system, and a mixture of water, 
oil and amphiphile is called a ternary system. In binary systems, amphiphiles 
self-assemble into spherical, cylindrical or bilayer structure in such a way that 
the hydrophobic tails are shielded from the hydrophilic solvent; in ternary sys- 
tems, amphiphiles self-assemble into monolayer structures in such a way that 
the hydrophobic tails and the hydrophilic heads face the hydrophobic and hy- 
drophilic solvents, respectively (see Fig. 2). Irrespective of these differences, 
amphiphiles in binary and ternary systems assemble into similar geometries, 
because in both cases one deals with extended sheet-like structures (the same 
geometries also occur in diblock copolymer systems). At room temperature, 
amphiphiles in mono- and bilayers usually form a two-dimensional fluid, that 
is the molecular constituents show only short-ranged and no long-ranged order 
in the plane of the interface. Therefore the sheet-like structures in amphiphilic 
systems are called fluid membranes (in fact membrane fluidity is essential for 
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Figure 3: Geometries of self-assembled interfaces: (a) micellar (spheres), (b) 
hexagonal (cylinders), (c) lamellar (planes), (d) cubic bicontinuous (cubic mini- 
mal surfaces) and (e) disordered bicontinuous (random surfaces). In amphiphilic 
systems, mono- or bilayers are draped onto the mathematical surfaces. Bicon- 
tinuous means that there is one surface which spans the whole sample, thereby 
separating it into two disjunct yet intertwined labyrinths (like the corresponding 
surface, each labyrinth is connected and spans the whole sample). Note that 
except for (e), all surfaces are ordered and have constant mean curvature. 

the functioning of the protein machinery carried by biomembranes) . 

Since the physics of amphiphilic systems is mostly determined by their inter- 
faces, different geometries have comparable free energies and small changes in 
external variables can induce phase transitions. In particularly, for amphiphilic 
systems different geometries are stable for different values of concentrations and 
temperature. In order to achieve systematic control of amphiphilic systems, a 
large effort has been invested into experimentally determining the phase dia- 
gram for many important amphiphilic systems. It has been found that despite 
their molecular diversity, the phase behavior of amphiphilic systems follows 
some general rules which result from their intrinsic tendency for self-assembly 
into sheet-like structures. For example, for increasing amphiphile concentra- 
tion one usually finds the generic phase sequence micellar disordered - micellar 
ordered - hexagonal - cubic bicontinuous - lamellar. Micellar means spherical 
geometry; micellar ordered corresponds to an ordered (usually cubic) array of 
spherical micelles. Hexagonal denotes a two-dimensional packing of cylindri- 
cal aggregates. Bicontinuous means that one surface partitions space into two 
separate labyrinths, each of which can be used to traverse space. Cubic bicon- 
tinuous corresponds to a space-filling arrangement of one interface, which folds 
onto itself in a cubic arrangement. Sometimes non-cubic spacegroups are found, 
but most ordered bicontinuous structures in equilibrium are cubic. In some sys- 
tems, bicontinuous yet disordered phases occur, which are called sponge phases 
and microemulsions for binary and ternary systems, respectively. Lamellar cor- 
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responds to a onc-dimcnsional stack of interfaces. Ordered arrangements like 
hexagonal or cubic symmetries can be probed by scattering techniques, and 
bicontinuity by diffusion experiments (e.g. using nuclear magnetic resonance). 
Fig. 3 tries to visualize the different interface geometries. It is the case of cubic 
bicontinuous phases which is the main focus of this paper. 

There exists a large body of literature on the occurrence of bicontinuous 
phases in amphiphilic systems (for reviews sec [27, 68, 103]). In Fig. 4 we show 
experimental phase diagrams for binary [72, 108] and ternary [59, 76] systems 
with the small surfactant (this is the molecule shown in Fig. la). For 

the binary system. Fig. 4a, at low temperature one sees the sequence hexagonal 
Hi - cubic bicontinuous Vi - lamellar La with increasing amphiphile concen- 
tration. For high temperature, the hexagonal and cubic bicontinuous phases 
disappear, the lamellar phase expands (its interfaces unbind), and for very 
small amphiphile concentration a sponge phase L3 occurs. For the ternary 
system, Fig. 4b, we recognize the same situation again close to the binary side 
watcr-amphiphilc. For equal amounts of water and oil and not too large concen- 
tration of amphiphile, that is in the lower middle of the Gibbs triangle, one sees 
the following phase behavior: at low temperature an emulsification failure oc- 
curs, that is the micellar disordered phase Li coexists with an oily excess phase. 
At high temperature, the emulsification failure disappears and is replaced by 
a microemulsion (the central, triangular one phase region). We can conclude 
that at low temperature, strongly curved structures prevail {Hi, Li), while at 
high temperature, those structures become stable which are locally flat {La, L3, 
microemulsion). As we will show below, this behavior can be understood using 
the concept of a temperature dependent spontaneous curvature: spontaneous 
curvature is finite at low temperature, but vanishes at high temperature. The 
concept of finite spontaneous curvature also offers a natural explanation for the 
emulsification failure at low temperature. As we will argue below, the cubic 
bicontinuous phase Vi is an intermediate structure between strongly curved and 
locally flat. It also follows from the phase behavior described that at low tem- 
perature, stable phases tend to be ordered {Hi, Vi), while disordered phases 
profit from increased temperature (L3, microemulsion). We will discuss below 
that the disordered phases indeed have a larger configurational entropy. In this 
respect, it is the lamellar phase La which acts as an intermediate structure. 

In Fig. 5 we show experimental phase diagrams for two different lipid-water 
mixtures. Fig. 5a is for water and 2:1 lauric acid / dilauroyl phosphatidylcholine 
[110] and Fig. 5b for water and monoolein [88] (these arc the molecules shown in 
Fig. Ib-d). Despite the molecular differences, the macroscopic phase behavior 
is surprisingly similar. At low temperature the membranes loose their fluidity. 
At intermediate temperatures, the lamellar phase is stable, and at high tem- 
peratures, it is replaced by the hexagonal phase. Whereas in the case of the 
surfactant Ci2E^ spontaneous curvature is finite at ambient temperatures and 
vanishes at high temperature, for the lipids spontaneous curvature vanishes at 
ambient temperatures and increases with temperature. Several cubic bicontin- 
uous phases are stable for intermediate temperatures, in the sequence lamellar 
- G - D - P (here G, D and P stand for cubic bicontinuous structures which 
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Figure 4: Experimentally determined phase behavior for the surfactant C12E5: 
binary system {H2O/C12E5) [108] (top) and ternary system {H2O/C12E5/C14) 
[76] (bottom). In the binary system, the hexagonal phase Hi, the cubic bicon- 
tinuous phase Vi, the lamellar phase La, the sponge phase L3 and the micellar 
phases Li and L2 are stable. In the ternary system, the same phases occur 
close to the binary side H2O/C12E5. At Ti = 5°T and T2 = 25°T, Li coexists 
with excess oil (emulsification failure). At T3 = 48° T (balanced temperature), 
a bicontinuous microemulsion is stable in the middle of the Gibbs triangle. 

are discussed in more detail below) with increasing water content. The last 
phase in this sequence undergoes an emulsification failure, that is it coexists 
with an excess water phase. We will show below that this phase behavior can 
be explained nicely if the lipid monolayers of the cubic bicontinuous phases are 
modeled as parallel surfaces to a cubic minimal midsurface. 

In this article, we focus on the interface description of amphiphilic systems, 
that is the free energy of the system is written as a functional of its inter- 
face configuration. The interface free energy is introduced in the next section 
(Sec. 2), and in the rest of this article we will specify this free energy expression 
for different instances of cubic bicontinuous phases. In each case, we will discuss 
the relevant geometric properties and show how they relate to the resulting free 
energy expressions and phase diagrams. Our starting point are triply periodic 
minimal surfaces (TPMS) in Sec. 3, a subject well-known from differential ge- 
ometry. These structures are expected to occur if the amphiphilic interface is 
symmetric in regard to its two sides and if temperature is not too high as to 
destroy the ordered state. TPMS are also the reference state for their parallel 
surfaces and CMC-companions, which are the adequate structural models for a 
detailed analysis of amphiphilic monolayers in cubic bicontinuous phases. In the 
case that the physical interface is not symmetric in regard to its two sides (that 
is if spontaneous curvature exists, like usually for an amphiphilic monolayer at 
a water-oil interface), the relevant mathematical representations are the con- 
stant mean curvature companions of the TPMS, that is triply periodic surfaces 
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Figure 5: Experimentally determined phase behavior for lipid-water mixtures: 
2:1 lauric acid / dilauroyl phosphatidylcholine [110] (top) and monoolein [88] 
(bottom). In both phase diagrams, one sees the sequence lamellar - cubic bi- 
continuous - excess water (emulsification failure) at intermediate temperatures 
(around 40°C) and with increasing water concentration. In (a), G, D and P are 
stable. In (b), G (Ia3d) and D (Pn3m) are stable. 

of constant mean curvature (CMC-surfaces) treated in Sec. 5. The case of lipid 
bilayers might be considered to be the composition of two such CMC-surfaces, 
but we will argue in Sec. 4 that in this case the relevant geometry is in fact 
the one of parallel surfaces to a TPMS. In all of these cases, we are interested 
in ordered structures, and therefore the main method will be minimization of 
the corresponding free energy functionals. This is different in Sec. 6, where we 
will discuss disordered bicontinuous phases, that is the sponge phases and mi- 
croemulsions, which often occur at higher temperatures due to entropic effects. 
These structures are modeled as random surfaces, thus the main model here 
will be Gibbs distributions, in particular the theory of Gaussian random fields 
and Monte Carlo simulations. 

2 Free energy functionals 

2.1 Interface models 

From the viewpoint of elasticity theory, amphiphilic interfaces with small cur- 
vatures can be considered to be thin elastic shells, which are known to have few 
fundamental modes of deformation: out-of-plane bending, in-plane compression 
and in-plane shearing [60]. However, since amphiphilic interfaces are fluid and 
nearly incompressible, in-plane strain is irrelevant and the most relevant defor- 
mation mode is bending. For small curvatures the free energy of an amphiphilic 
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interface is a function only of its geometry [8, 47]: 

F = J dA{a + 2K{H-cof + RK} . (1) 

Here dA denotes the differential area element and H and K mean and Gaussian 
curvature, respectively. The latter two follow from the two principal curvatures 
fci and ^2 as _ff = (fci + fc2)/2 and K = fcifc2. The three material parameters 
introduced in Eq. (1) define the energy scales of the corresponding changes: 
a is surface tension and corresponds to changes in surface area, k is bending 
rigidity and corresponds to cylindrical bending, and R is saddle-splay modulus 
(or Gaussian bending rigidity) and corresponds to changes in topology due to 
the Gauss-Bonnet-theorem for closed surfaces, / dAK = 2-kx, where x denotes 
Eulcr characteristic. The spontaneous curvature defines the reference point for 
bending deformations. It has to vanish for symmetric amphiphilic sheets (like 
lipid bilayers), but in general is finite for non-symmetric ones (like monolayers). 
For monolayers, it is usually a linear function of temperature, cq ~ (T — T;,). 
Therefore spontaneous curvature vanishes at the balanced temperature T^, at 
which solvent properties make the monolayer symmetric in regard to bending. 

For amphiphilic interfaces, surface area is proportional to the number of 
amphiphilic molecules, thus a can also be interpreted as a chemical potential 
for amphiphiles. The bending rigidity k has to be positive, otherwise the system 
would become instable to spontaneous convolutions. For surfactant and lipid 
systems, its values are of the order of 1 and 20 ksT, respectively. Surface 
tension u can assume negative values (favorable chemical potential for the influx 
of amphiphiles). as long as bending rigidity n exists and has a positive value, 
in order to prevent an instability towards proliferation of interfacial area. The 
value of the saddle-splay modulus R is often debated, but usually it is assumed 
to have a small negative value. This assiimption is validated by the following 
argument [48]: for a = and cq = 0, we can rewrite Eq. (1) as 

F = j dAi^K+{kt+k2f + \n-{ki-k2f^ (2) 

with 

= 2' '^- = ~2' ^' 
Prom this we conclude that the Gaussian bending rigidity has to satisfy k+ > 

and K- > 0, that is —2k < k < 0, otherwise the system would become unstable. 
For R < —2k (k+ < 0), we would get ki = k2 ^ oo, that is many small droplets, 
and for K > (k_ < 0), we would get fci = — fc2 oo, that is a minimal surface 
with very small lattice constant. 

From the mathematical point of view, it is interesting that for a = and 
Co = 0, the bending energy from Eq. (1) is not only invariant under rescaling with 
length, but also invariant under conformal transformations in general. This has 
intriguing consequences for vesicles [104] (vesicles are depicted in Fig. 2, but are 
not subject of this article) and phase transitions in systems without spontaneous 
curvature (as will be discussed in Sec. 6). 
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Since the energy scales involved are of the order of ksT, thermal noise is suf- 
ficient to induce shape changes. Due to thermal fluctuations on smaller length 
scales, the effective values for the material parameters are changed (renormal- 
izcd) at larger length scales [49, 79]. It has been shown in the framework of 
renormalization group theory that logarithmic corrections arise due the two- 
dimensional nature of the amphiphilic interfaces, so that [79, 17, 7] 

'SksT I 

K^(/) = «-__ln- , (4) 
kr{1) = k+ — — In - (5) 

OTT 

where I is the length scale on which the system is analyzed and 5 is a microscopic 
cutoff on the scale of the membrane thickness. Thus bending rigidity decreases, 
while Gaussian bending rigidity increases with increasing length scale I. 

For the following, it is useful to introduce two dimensionless quantities which 
correspond to the two curvatures of a two-dimensional surface. For a structure 
with surface area A and volume V, the length scale V/A can be used to rescale 
its curvatures. In order to be able to convert easily from local to global quanti- 
ties, we consider a triply periodic CMC-surface with Eulcr characteristic x and 
integrated mean curvature Hi = J dAH per conventional (that is simple cubic) 
unit cell: 

These expressions motivate the definition of the topology index T and the cur- 
vature index A: 

r=(^y A = ^ (7) 

where A* = AjV'^l'^ and H* = Hi/V^^'^ arc the scaled surface area and the 
scaled integral mean curvature per conventional unit cell, respectively. Note 
that these definitions can be applied to any triply periodic surface; here we used 
CAf C-surfaccs only for an heuristic motivation of the quantities defined. Both 
the topology index F and the curvature index A do not depend on scaling and 
choice of unit cell. They are universal geometrical quantities which characterize 
a surface in three-dimensional space. The topology index describes its porosity 
(the larger its value, the less holes the structure has) and its specific area content 
(the larger its value, the more inner surface the structure contains), and the 
curvature index describes how strongly the structure is curved (irrespective of 
the actual lattice constant). It is interesting to note that the two indices defined 
here correspond to the two isoperimetric ratios known from integral geometry. 
For a TPMS, the curvature index A vanishes and the topology index T is its 
most important geometric characteristic. 
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2.2 Ginzburg-Landau models 

A different but equally powerful approach to ampliiphilic interfaces is the iso- 
surface construction (also known as phase field method), which derives a two- 
dimensional surface from a three-dimensional scalar field 3>(r). For ternary 
amphiphilic systems, $(r) can be interpreted as the local concentration differ- 
ence between water ($ = 1) and oil ($ = —1). The position of the amphiphilic 
monolayer can be identified with the isosurface $(r) = 0. For binary amphiphilic 
systems, $(r) can be interpreted as the local concentration difference between 
water on different sides of a bilayer, and the isosurface $(r) = marks the 
position of the bilayer mid-surface. In both cases, an energy functional can be 
defined in the spirit of a Ginzburg-Landau theory for the scalar field $(r). A 
reasonable choice is the $^-model introduced by Gompper and Schick [39, 40]: 

= Jdr {(A$)2 + gmW^f + /($)} (8) 

with the following choice for / and g: 

/($) = ($ + - 1)2($2 + /q) , ^($) = 50 + 92$' . (9) 

It has been shown that this model is similar to an interface description as given 
in Eq. (1), and prescriptions have been given how to calculate the parameters of 
the interface Hamiltonian from the given Ginzburg-Landau theory [41] . Due to 
the invariance under $ — > — spontaneous curvature cq vanishes in this model. 
For given model parameters {go, 92, /o)j this energy functional can be minimized 
for its spatial degrees of freedom, and a phase diagram can be calculated as a 
function of model parameters by idcuitifying the absolute minimum at every 
point. A reasonable choice for amphiphilic systems is go < 0, fo close to and 
52 not too large. Then the formation of interfaces is favored and a lamellar phase 
becomes stable. Calculation of the corresponding parameters of the interface 
model then yields a < 0, k > and R <0. Note that a < favors the formation 
of interfaces, and that k < favors the lamellar phase. 

3 Triply periodic minimal surfaces (TPMS) 

Minimal surfaces and surfaces of constant mean curvature (see Sec. 5) in general 

have attracted a lot of attention both in mathematics and in physics, partially 
due to their intrinsic beauty, which might be considered to follow from the fact 
that they are solutions to variational problems [51]. It is important to note 
that the different fields are interested in different aspects of the same object: 
while for mathematicians the most central question is the existence proof for 
a minimal surface of interest, physicists are more interested in its representa- 
tion, which can be used to derive physical properties of corresponding material 
systems. As Karcher and Polthier remark, outside mathematics only pictured 
minimal surfaces have been accepted as existent [56]. In fact this statement 
can also be reversed: not every structure which one can picture is necessarily 
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a minimal surface. In this article we consider only established minimal sur- 
faces. As we are motivated by physical considerations, we are interested only in 
embedded surfaces which do not self- intersect. In regard to mathematics, our 
main emphasis here will be on representations, since the physical properties of 
amphiphilic systems follow from their spatial structure. 

Minimal surfaces are surfaces with H = everywhere. If one varies a sur- 
face with a normal displacement 5(j}{u^v) (where u and v arc the internal co- 
ordinates of the surface), it can be shown that the change in surface area is 
= 25 J dA(p{u,v)H{u,v) + 0{6'^) [54], that is a minimal surface is a sta- 
tionary surface for variations of surface area. Therefore surfaces under surface 
tension (like soap films), which try to minimize surface area, form minimal sur- 
faces. Here we consider surfaces which are dominated by bending rigidity rather 
than by surface tension. However, for vanishing spontaneous curvature the re- 
sulting structures also correspond to minimal surfaces. Eq. (1) explains why: 
these surfaces have to minimize / dAH'^ (they are so-called Willmore surfaces), 
and minimal surfaces arc a special case of Willmore surfaces, since iJ = is a 
trivial minimization of this functional. 

If mean curvature H = {ki+k2)/2 = 0, then Gaussian curvature K = k\k2 
■ < 0. Thus minimal surfaces are everywhere either saddle- like or flat, and 
nowhere convex. In fact it can be shown that the flat points with ki = k2 = 
are isolated. Moreover it follows from K < that minimal surfaces without 
boundaries cannot be compact [54]. Until 1983, only two embedded minimal 
surfaces were known which are non-periodic; the plane and the catenoid. Then 
a new surface of this type was found, the Costa-surfaces, which from the dis- 
tance looks like the imion of a plane and a catenoid [56]. Today some more 
surfaces of this type are known, but the majority of all embedded minimal sur- 
faces without boundaries are in fact periodic. There is only one simply-periodic 
minimal surface, the helicoid, and there are some doubly-periodic surfaces, with 
the Scherk-surface being most prominent. The majority of all known periodic 
minimal surfaces is triply-periodic. Here we will focus on cubic minimal surfaces 
and their physical realizations, which are cubic bicontinuous structures. 

The best known triply periodic minimal surface (TPMS) is the P-surface 
found by Schwarz in 1867 and depicted in Fig. 3d. Schwarz and his students 
found 5 TPMS, including the cubic cases P, C(P) and D. Until 1970 no more 
examples were found, then Schoen described 13 more [96], including the cu- 
bic cases G, F-RD, I-WP, 0,C-TO and C(D). He proved the existence of G 
by providing an explicit (Weierstrass) representation, and the existence of the 
others was proven in 1989 by Karcher [55]. More TPMS have been discovered 
by Fischer and Koch [23, 57] and others, and today many more can be gener- 
ated by making controlled modifications to computer models of known TPMS, 
a method which was pioneered by Polthier and Karcher [56]. However, these 
surfaces tend to be rather complicated, and we will show below that only the 
simple (as quantified by the topology index) TPMS are relevant for amphiphilic 
systems. In Fig. 6 we show visualizations of 10 different TPMS with cubic sym- 
metry. There representations have been obtained from the $^-Ginzburg-Landau 
theory as will be explained below. For some of these structures, we also show 
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Figure 6: Visualizations of 10 different cubic minimal surfaces. Every triply 
periodic surface divides space into two labyrinths, which in some cases are rep- 
resented here by skeletal graphs. D and G are the basis for the double-diamond 
structure and the gyroid, respectively, which often occur in physical systems. 

line-like representations of the two labyrinths defined by the surface {skeletal 
graphs). For example, in the case of the Schwarz P-surface, the two skeletal 
graphs are the edges of the surrounding cube and the three lines through the 
origin. 

Different methods can be used to obtain representations of TPMS. Until 
recently, the main method were the Weierstrass representation formulae. For the 
cubic TPMS, they are known for P, D and G [25, 26] as well as for I-WP [62, 12]. 
For each of these TPMS, a fundamental domain can be identified, so that the 
rest of the surface follows by replicating it with the appropriate space group 
symmetries (/m3m, Pn3m, laid and ImSm, respectively). The Weierstrass 
representation is a conformal mapping of certain complicated regions within 
the complex plane onto the fundamental domain: 



where (w, v) are the internal (and conformal) coordinates of the minimal surface. 
The Weierstrass mapping can be understood as the inversion of the following 
composition: first the surface is mapped onto the unit sphere via its normal, 
and then the unit sphere is mapped onto the complex plane by stereographic 
projection. The geometrical properties of a surface follow from the Weierstrass 
representation as 



dAiz) = \R{z)f (1 + |z|2)2 dudv, H{z) = 0, K{z) = , ..^.i ^ , ,2U (H) 



with z = u + iv. Obviously the (isolated) poles of R{z) correspond to the flat 
points {K = 0) of the minimal surface. Only few choices of R{z) yield embedded 
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X 


A* 


F 


G 


-8 


3(1 + k:()/2ki = 3.091444 


0.766668 


D 


-1() 


3/A:i = 3.837785 


0.749844 


I-WP 


-12 


2V3 = 3.464102 


0.742515 


P 


-4 


Sfci = 2.345103 


0.716346 


C(P) 


-16 


2,/k2 = 3.510478 


0.655993 



Tabic 1: Eulcr characteristic scaled surface area A* in the conventional unit 
cell and topology index F = (A*^/27r|x|)^/^ for those TPMS, for which exact re- 
sults are known from Weierstrass representations. Here ki = K{l/2)/K{\/i/2) 
where K{k) is the complete elliptic integral of the first kind, and ^2 = 
K{l/^/2,)/ K[^2/2,). Note that often a unit cell is chosen for D which is the 
eighth part of the one chosen here; then one has x = — 2 and A* = 1.918893. 

minimal surfaces. The ones for D and P have been known since the 19th century 

from the work of Schwarz: for D it is R{z) = {z^ — 14z^ + 1)~ 2 . P follows simply 
by the Bonnet transformation R{z) e^^R{z) with 6 = 90°. Eq. (11) implies 
that P and D have the same metric and the same distribution of Gaussian 
curvature. However, since they map differently into embedding space, they 
have different space groups and lattice constants. The gyroid G was discovered 
in 1970 by Schoen [96] as another Bonnet transformation of D, with = 38.015°. 
The Weierstrass representation for I-WP was found only recently [62, 12]. If one 
of its poles is chosen to be at infinity, one has R{z) = {z{z^ + l))-3. In Tab. 1 
we give exact results for geometrical properties which have been derived from 
Weierstrass representations (although for C(P) no Weierstrass representation is 
known, these values can be derived from its complementary relationship to P). 
Note that the topology index F establishes the hierarchy G - D - P within this 
Bonnet family. This sequence corresponds to the connectivity of the labyrinths 
defined by the TPMS: G has 3-fold, D 4- fold and P 6- fold coordination, since 
higher coordinated structures have more holes. 

It is interesting to note that the identification of these structures in material 
systems is a difficult enterprise, which depends on the availability of suitable rep- 
resentations. Physical representations of TPMS have been built since the 19th 
century, using soap films draped onto wire skeletons (without the wire skeletons 
surface tension would shrink these structures into collapse). In 1967 Luzzati and 
Spegt noted that certain phases in lipid systems are cubic bicontinuous [67], but 
it was only in 1976 that Sciven suggested the relevance of TPMS as structural 
models for cubic bicontinuous phases [102]. The Luzzati-Spegt structure was 
later identified with the gyroid structure G, and the diamond structure D was 
identified both in lipid [65] and diblock copolymer systems [114]. For diblock 
copolymer systems, the gyroid structure G was identified in 1994 [45], and little 
later it was noted that earlier identifications of cubic bicontinuous phases in di- 
block copolymer systems often mistook the gyroid structure G for the diamond 
structure D [46] (in fact, for diblock copolymer systems one should rather use 
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the terms double diamond and double gyroid, since in this case two interfaces are 
arranged around the corresponding TPMS). Today, the structures G, D, P and 
I-WP seem to be non- ambiguously identified in amphiphilic system, with exper- 
imental evidence based mainly on small angle scattering experiments, electron 
transmission microscopy, and swelling and diffusion experiments [27, 68, 103]. 

Apart from the cases G, D, P and I-WP given above, no more Weierstrass 
representations arc known for cubic TPMS, so for all other cases numerical meth- 
ods have to be used. For example, the Surface Evolver is a software package 
written by Brakke which allows to minimize triangulated surfaces for different 
energy functionals (compare Sec. 8). Here we discuss the method of constructing 
the isosurfaces of a scalar field $(r), as introduced in Sec. 2.2. For our purpose, 
the usefulness of this model lies in its rugged energy landscape, which means 
that many more local minima exist than the absolute minima corresponding to 
the lamellar phase. In physical terms, these additional minima correspond to 
modulated phases which are metastable. If started with suitable initial con- 
ditions, the minimization procedure therefore yields representations which can 
be used to characterize the structural properties of these phases. Therefore 
this model has been used repeatedly in order to investigate bicontinuous cubic 
phases [41, 42, 43, 97]. In particularly, it has been found that the resulting 
representations are very close to TPMS [42, 43]. This finding can be explained 
as follows [97]: for a triply-periodic cubic structure, the free energy per unit 
volume / = F/V follows from the interface description of Eq. (1) as 



/ = i (aA*) + ^(2kJ dA + 2R7TX ) (12) 



where a is the lattice constant of the conventional unit cell. Both terms in 
brackets are scale invariant, that is they do not depend on the lattice constant 
a. Since for the Ginzburg-Landau model at hand the first and second term 
in brackets is negative and positive, respectively (compare Sec. 2.2), a balance 
exists between the negative surface tension term, which favors small values 
of a, and the positive curvature contributions, which favors large a. If we 
make the assumption that the minimization problem can be decomposed in two 
independent minimizations, one for lattice constant a and one for the shape, we 
can first minimize / for a while considering A* and / dA to be constant: 

) ' -/"---(^^j [2RTTx + 2KjdAH'^ 

(13) 

Similar mechanisms are always at work in amphiphilic systems and explain why 
in contrast to surface tension, bending rigidity does not necessitate the presence 
of a scaffold (like a wire skeleton for soap films) to prevent collapse. In a second 
step, we now minimize fmin for shape. If we assume A* to be constant, we have 
reduced our problem to the Willmore problem, and this is why minimal surfaces 
with H = appear in the framework of the $^-theory. The final free-energy 
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density can be written as 

fmin = - (27) (\^) r 

where F is the topology index defined in Eq. (7). We conclude that the stability 

of the different TPMS is governed by their topology indices and that the most 
favorable bicontinuous cubic phase is the gyroid G since it has the largest value 
for r (compare Tab. 1). Although this conclusion is based on some assumptions, 
it is corroborated by a detailed numerical analysis [97] . 

In order to obtain representations for a large list of TPMS, we used the 
Fourier approach and the theories of black and white space groups to obtain 
TPMS as local minima of the $^-theory [97]. Black and white space groups 
are also known as magnetic or Shubnikov space groups and lead to a crystal- 
lographic classification of TPMS [23, 24]. Implementation of black and white 
symmetries leads to a considerable reduction in the degrees of freedoms of the 
Fourier series. Therefore this method is computationally cheap and its results 
are easy to document and to reuse. A triply periodic surface partitions space 
into two labyrinths, which can be considered to be colored black and white. In 
the framework of the Ginzburg-Landau theory, black and white correspond to 
$ > and $ < 0, respectively. The surface is called balanced if there exists an 
Euclidean transformation a which maps the white labyrinth onto the black one 
and vice versa; otherwise it is called non-balanced. Examples for non-balanced 
surfaces are I-WP and F-RD, and the Fourier approach for the corresponding 
structure follows the usual rules for the respective space group [106]. However, 
if the surface is balanced, the structure is characterized by two space groups: 
if the colored structure has space group H, the uncolored structure has space 
group Q = Ti. (E) {1, a}. Ti contains all symmetry operations of Q that do not 
interchange the two labyrinths. It is a subgroup of Q of index 2, that is the 
quotient group is isomorphic to the cycHc group of order 2, Q /H = Z2 = {1, a}. 
For example, for the balanced cubic surface defined by cos x + cos y + cos z = 0, 
the operation interchanging the two labyrinth is the translation by half a body 
diagonal, and we have 7i = Pm3m and G = ImSm. 

Since 7i has index 2 in Q, it follows from the theorem of Hermann [105] 
that 7i has cither the same point group or the same Bravais lattice as Q. If Ti 
and Q have the same point group, their Bravais lattices have to be different. 
Thus the Euclidean operation a. has to be a translation which extends one cubic 
Bravais-lattice into another. For the cubic system, there are only three different 
Bravais lattices (simple cubic P, body-centered cubic I and face-centered cubic 
F), and only two possibilities to extend one into the other by a translation to,: 
for toe = + y + z)/2 a P-latticc becomes a I-lattice, and for t^ = ax/2 
a F-lattice becomes a P-lattice. The condition <I>(r + ta) = ^ *I'(r) leads to 
reflection conditions for the reciprocal vectors. For Miller indices {h,k,l), one 
finds /i -I- fc -f / = 2n + 1 for P ^ /, and h, A;, / = 2n + 1 for F ^ P. Note that 
these reflection conditions are similar to the well-known ones h + k + I = 2n for 
I and h + k, h + 1, k + 1 = 2n for F . The P-surface has space group Pm3m, but 
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n 


Q 




nodal approximations 


G 


/4i32 


Ia3d 


432 m3m 


sin(.x) cos(y) + sin(y) cos(z) + sin(;z) cos(a;) 


D 




Pn3m 


P 


cos(a; — y) cos(z) + sin(a:; + y) sin(z) 


P 


Pm?>m 


ImSm 


P^ I 


cos(a;) + cos(y) + cos(2;) 


C(P) 


Pm2>m 


Im3m 


P^ I 


cos(a;) + cos(y) + cos(^) + 3 cos(x) cos(t/) cos(2;) 



Table 2: Group - subgroup pairs G — H, the relation H ^ G between them 
and nodal approximations for some balanced cubic minimal surfaces. Q and 
Ti. differ cither in Bravais lattice or in point group, but not in both. Nodal 
approximations only consider the space group information given hy G — Ti-. For 
the simple cases, the first mode of the corresponding Fourier series is sufficient 
to obtain the correct topology. 





G 


nodal a.pproxiiui 1 1 i( )iis 


1-WP 


Im3m 


2[cos(a;) cos(y) + cos(y) cos(z) + cos(2:) cos(x)] 
— [cos(2a;) + cos(2y) + cos(22;)] 


F-RD 


FmSm 


4 cos(a;) cos(j/) cos{z) 

— [cos(2a;) cos(2?/) + cos(2y) cos(22;) + cos(2^;) cos(2a;)] 



Table 3: Nodal approximations for some non-balanced cubic minimal surfaces, 
for which H = g. In these cases, more than one mode is needed. 

a black and white symmetry with supergroup ImSm, therefore it is an example 
for the case P ^ I. Thus we have h + k + I = 2n + 1 and the first Fourier 
mode is (1, 0, 0). Note that the double P structure, which one would expect for 
a diblock copolymer system, has no black and white symmetry, but space group 
Im,3m,. thus one has h + k + I = 2n and the first mode is (1, 1, 0). 

If Ti. and G have the same Bravais lattice, their point groups have to be 
different. Thus the Euclidean operation a has to be a point group operation 
which extends one cubic point group into another one. There are five cubic 
point groups and six ways to extend one of them into another by some Pa ■ The 
condition <S>{Par) = — ^(r) loads to more complicated rules than in the case of 
identical points groups [97]. However, there are few relevant examples from this 
class, the most interesting one being the gyroid G, where Pa is the inversion and 
the resulting rules are rather simple again (the even part of the Fourier series 
vanishes). 

Fischer and Koch have completely enumerated all 34 cubic group-subgroup 
pairs G — "H with index 2 which are compatible with cubic balanced TPMS 
[23]. Although there is no way to completely enumerate all TPMS belonging 
to a given pair C/ — ?i in general, this is possible for a certain subset, that is 
for all cubic balanced TPMS which contain straight lines which form a three- 
dimensional network [23]. This list reads P, C(P), D, C(D), S and C(Y). In 
our work, we implemented Fourier series for these structures as well as for the 
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Figure 7: Distribution f{K) of Gaussian curvature K over the cubic minimal 
surfaces investigated in the framework of the $^-model. (a) For the structure 
G, D, P and I-WP, the data is obtained from their exact Weierstrass representa- 
tions; the agreement with the numerical results from the $^-model is good (not 
shown), (b) For the structures F-RD, S and C(P), no Weierstrass representation 
is known, and the data is obtained from the <f>^-model. 



gyroid G (which is balanced, but does not contain any straight lines) and the 
non-balanced structures I-WP and F-RD. In some cases like G, D and P, the first 
mode of the Fourier ansatz implementing the correct black and white symmetry 
is already sufficient to obtain a representation which is topologically correct 
{nodal approximation). In all other cases investigated, the addition of one more 
mode (with the relative weight fixed by visual inspection) is sufficient. Nodal 
approximations for bicontinuous structures have been discussed first by von 
Schnering and Nesper [115]. In Tab. 2 and Tab. 3 we give nodal approximations 
and the corresponding space group information for some of the balanced and 
non-balanced structures investigated. 

Implementing the complete Fourier series and numerically minimizing the 
$^-functional from Eq. (8) with nodal approximations as initial conditions, we 
arrived at improved nodal approximations which were tabulated in Ref. [97] with 
up to six Fourier modes. The isosurface construction is easily implemented us- 
ing the marching cube algorithm. The resulting triangulations can then be used 
to investigate geometrical properties of these surfaces. In particular, widely 
used mathematics programs (like Mathematica, Maple or Matlab) can be eas- 
ily used to obtain useful representations of TPMS from (improved) nodal ap- 
proximations. Our numerical work shows that the six modes of the improved 
nodal approximations are sufficient to decrease the deviation of the curvature 
properties from the real TPMS by one order of magnitude compared to the 
nodal approximations. Also we measured for the first time the distribution 
f{K) = J dA{u,v) 5{K — K{u,v)) of Gaussian curvature K over the surfaces. 
For this purpose, we used up to 100 Fourier modes. The distributions f(K) 
for different structures are plotted in Fig. 7. In the cases for which Weierstrass 
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G, D, P 


I-WP 


S 


F-RD 


C(P) 


A 


0.218702 


0.482666 


0.586079 


0.649801 


0.842022 



Table 4: Variance A of the distributions of Gaussian curvature f{K). The 
vahics for A arc the same for G, D and P due to the existence of a Bonnet- 
transformation between them. 



representations are known, the same data has been derived from the exact rep- 
resentations, and the resulting agreement was very good. In general, we found 
that the surfaces C(D), C(P), F-RD, S and C(Y), for which no Weierstrass rep- 
resentations are known, are much more compHcated than the cases G, D, P and 
TWP. This is evidenced by larger values of the topology index F, multi-modal 
distributions f{K) of Gaussian curvature K, and the larger widths of these dis- 
tributions. The latter property can be quantified by defining a dimensionless 
variance 

where (...) means area average. In Tab. 4 we give the variance A of the dif- 
ferent distributions. Its value is the same for G, D and P due to the existence 
of a Bonnet-transformation between them. This means that G, D and P have 
the narrowest distribution of K (they are most uniformly curved), and all other 
structures have a considerably wider one, with I-WP being the next best struc- 
ture. 



A = A^11^^ (15) 



4 Parallel surfaces 

We now consider the case of cubic bicontinuous phases in lipid-water mixtures. 
The two examples for the experimental phase behavior of such systems given 
in Fig. 5 show obvious similarities, and we will show now that a theoretical 
description can nicely explain these [99, 100]. It has been shown by a thorough 
analysis of electron density maps derived from X-ray data that the mid-surfaces 
of the cubic bicontinuous structures in these systems are very close to minimal 
surfaces [68]. Indeed if one considers the lipid bilayer as one entity, it has no 
spontaneous curvature by symmetry and one expects a minimal surface shape 
for the midsurface. Each of the two monolayers of the bilayer has unequal sides 
and therefore finite spontaneous curvature. For lipids monolayers, spontaneous 
curvature increases towards the water side as a linear function of temperature. 
From the interface point of view, one expects each monolayer to form a CMC- 
surface. However, then the distance to the minimal mid-surface would vary 
with position and the amphiphilic tails would have to stretch in order to fill 
the internal space of the membrane [3, 10]. It can be shown in the framework 
of a simple microscopic model that the relative importance of stretching to 
bending contributions to the free energy of the bilayer scales as {a/S)"^, where a 
is lattice constant and 5 is tail length [100]. Therefore stretching is prohibitively 
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expensive and the two monolayers do not form CMC-surfaces as expected from 
the curvature energy of Eq. (1), but rather paraUel surfaces to the minimal 
mid-surface. The free energy of the cubic bicontinuous phases then follows by 
specifying Eq. (1) with spontaneous curvature for two monolayers which arc 
parallel surfaces to a given TPMS. Since the parallel surface geometry does 
not allow to completely relax the bending energy, the overall structure is called 
frustrated [3, 10, 20]. 

In principle, the analysis in the framework of the parallel surface model is 
rather simple, since there exist exact formulae which express the geometrical 
properties of the parallel surface as a function of the geometrical properties of 
the minimal surface: 

dA^ = Mil+ K6^) , = ^ , = ^ (16) 

where S is the distance between the two surfaces (that is amphiphilic chain 
length). Note that these formulae are an extension of Steiner's theorem from 
integral geometry to non-convex bodies, which is valid as long as the distance 5 is 
smaller than the smallest radius of curvature of the surface (otherwise the surface 
will self-intersect). Since mean curvature H vanishes on the reference surface, 
the bending energy now becomes a function only of its Gaussian curvature K. 
However, as we have seen above, K is distributed over the TPMS in a non-trivial 
way, which makes the detailed analysis rather complicated. Nevertheless, one 
central conclusion can be already made at this point: using Eq. (16) in Eq. (1) 
leads to an effective curvature energy for the lipid bilayer (to second order in 5) 
[85] 

Fbi = j {4:cl6^K + (2k + 8co6k + Ac16^k)K + AkS^K^} . (17) 

Thus we see that although effective spontaneous curvature vanishes due to the 
bilayer symmetry, the effective saddle-splay modulus, Kbi = 2R. + 8cq6k + 0{S'^), 
is corrected to higher positive values due to the presence of the monolayer spon- 
taneous curvature cq. We conclude that as long as cqS > —k/Ak, the bicon- 
tinuous cubic phases are favored over the lamellar phase, since the preferred 
curvature of the monolayers translates into a topological advantage of saddle- 
type bilayer structures. We also see that the correction term in Rbi is linear 
in Co and therefore in temperature T. Therefore bicontinuous phases become 
more favorable with increasing temperature in general (compare also Sec. 6 on 
disordered bicontinuous structures). 

For a detailed analysis of cubic bicontinuous phases made from bilayers, we 
first note that the volume fraction of the lipid tails (the hydrocarbon volume 
fraction, which for simplicity we identify with the lipid volume fraction, since 
the lipid heads are rather small) can be calculated as 

where again we have used the Gauss-Bonnet theorem. Eq. (18) can be inverted 
to give a/ 5, the lattice constant a in units of the chain length S, as a function 
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Figure 8: For cubic bicontinuous phases in lipid-water mixtures, the hpid mono- 
layers can be modeled as parallel surfaces to a minimal mid-surface, (a) Free 
energy densities / as a function of lipid volume fraction v for several cubic bi- 
continuous phases and the lamellar phase, (b) Theoretical phase diagram as a 
function of water volume fraction and (dimensionless) spontaneous curvature. 



of hydrocarbon volume v. For small v we find 

and one can check numerically that a/6 = 2A* /v is an excellent approximation 
for V < 0.8. For larger values of v, the surfaces begin to self-intersect and 
our model becomes unphysical. Combining Eq. (1) and Eq. (16) yields for the 
bending energy per unit volume of the two monolayers (we use a factor S/AkCq 
to write this expression dimensionless and a factor 6 to write cq dimensionless) 




where K* = Ka^ is scaled Gaussian curvature, 5/a has been replaced by v/2A* ^ 
S(i4r*) has been defined as K* A* /%ttx (this can be considered to be a local 
analogue of F^), and r as — k/2k (0 < r < 1 due to the restrictions on R). The 
free-energy density / now is a function of lipid volume fraction w, spontaneous 
curvature cq, the ratio r of the two bending constants, and the distribution 
f{K) of Gaussian curvature K . The (numerical) evaluation of this expression 
is only possible with the knowledge of the f{K) for all TPMS of interest, which 
have been numerically obtained in Ref. [97] as described above. 

In the limit of a planar mid-surface (lamellar phase), Eq. (20) simplifies to 
f = V - the larger the lipid volume fraction, the more frustrated bending energy 
per volume accumulates. For a full analysis, one also has to include the effect 
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of thermal fluctuations. For the lamellar phase, they lead to a sterie repulsion 
between the interfaces, which lead to an additional term ~ — v)"^ for the 

lamellar phase. For the cubic bicontinuous phases, steric repulsion is irrelevant 
since the lateral restriction on the scale of a lattice constant leads only to small 
perpendicular excursions. However, here the renormalization of the saddle splay 
modulus K becomes relevant, and the formula given in Eq. (5) has to be incor- 
porated. The renormalization of bending rigidity k and spontaneous curvature 
Co is irrelevant here, since thermal fluctuations occur mainly on the level of the 
lipid bilayer, for which mean curvature H vanishes. For the lamellar phase, 
Gaussian curvature K vanishes as well, and the renormalization of all mate- 
rial parameters is irrelevant. Putting everything together, we can numerically 
calculate phase diagrams from the free energy densities of the different phases 
by using the Maxwell construction (construction of convex hull). In Fig. 8 we 
show both the free energy densities as a function of lipid volume fraction v (for 
fixed spontaneous curvature) and the theoretical phase diagram as a function 
of water volume fraction and spontaneous curvature. The main results are in 
excellent agreement with the experimental phase diagram shown in Fig. 5: from 
the many TPMS considered, only G, D and P are stable, their regions of sta- 
bility have the shape of shifted parabolae, and they occur in the sequence L - 
G - D - P - emulsification failure. These results can be understood as follows: 
G, D and P can achieve the least frustration since they have the narrowest dis- 
tribution of Gaussian curvature as measured by the variance A given in Tab. 4. 
This prediction has been stated before by Helfrich and Rennschuh [50], but at 
that time hardly any data was known to support it. The lipid volume fraction 
at which they achieve this can be estimated by setting the dimensionless mean 
curvature averaged over the parallel surface 

JdAS -4-{v/Tr ^^^^ 

equal to the spontaneous curvature cq. Therefore cq as a function of v essentially 
scales as ~ (1 — pw)"^ f^"^, where pw = 1 — v is the water volume fraction. This 
explains the characteristic shape of the bicontinuous stability regions in both the 
theoretical and experimental phase diagrams. Note that at high temperature, 
that is large spontaneous curvature, eventually the hexagonal phase will become 
stable, which cannot be treated in the framework presented here. The optimal 
lipid volume fraction v follows from Eq. (21) as 



4co 
1 + co 



r . (22) 



Therefore the structures G, D and P become stable in the sequence of their 
geometry index F, that is as G - D - P. If the water content corresponds to a 
average curvature of P which is larger than the optimal (spontaneous) curvature, 
some of the water is simply expelled from the structure, in order to keep the 
optimal curvature (emulsification failure). Finally it should be noted that the 
Bonnet-transformation connecting G, D and P causes the three structures to be 
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Figure 9: Geometrical data for triply periodic surfaces of constant mean curva- 
ture: (a) scaled surface area A* and (b) scaled mean curvature H* as a function 
of volume fraction v. The two branches are symmetrical for G, D, P and C(P) 
since their minimal surface members are balanced. 

stable along a triple line (see inset of Fig. 8). This means that the stability of D 
and P is very delicate and can easily be destroyed by additional physical effects. 
Therefore we conclude that in the framework of our interfacial approach, from 
all TPMS considered the gyroid G is the only phase which has a robust stability 
in lipid-water mixtures, since among the phases with favorable distribution of 
Gaussian curvature, its geometrical properties are closest to the ones of the 
lamellar phase. 

5 Surfaces of constant mean curvature (CMC) 

Surfaces of constant mean curvature (CMC-surfaces) have H — const every- 
where, so minimal surfaces are a special case of CMC-surfaces. In contrast to 
minimal surfaces, CMC-surfaces with finite H can be compact, but the only 
CMC-surface which is compact and embedded is the sphere (for a long time, 
the sphere was believed to be the only compact CMC-surface, but in 1986 Wente 
found the first CMC-torus) . Non-compact embedded CMC-surfaces are cylinder 
and unduloid, and all other known CMC-surfaces are doubly or triply periodic. 

CMC-surfaces are solutions to the variational problem of minimal surface 
area under a volume constraint. This can be shown as follows: we first introduce 
a Lagrange parameter for volume, that is pressure p. The corresponding energy 
is —pV. For normal variations S(p{u, v), we have AF — 2Sa J dA(t){u, v)H{u, v) — 
p5 J dA(j){u, v) + 0{S^). If the surface is required to be stationary in regard to 
variations in S, we obtain the Laplace equation H = p/2a and H is constant 
over the whole surface. This explains why soap bubbles and liquid droplets are 
spheres (respectively spherical caps when bound by a surface). In amphiphilic 
systems, CMC-surface arise in the presence of spontaneous curvature: like min- 
imal surfaces minimize the Willmore functional / dAH^ , CMC-surfaces with 
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G 


D 


I-WP 


P 


F-RD 


C(P) 


inc^asurcid 


0.2191 


0.1411 


0.1385 


0.2117 


0.0665 


0.0466 


api)roximutioii 


0.1901 


0.14G.J 


0.1592 


0.2188 


0.0906 


0.1226 



Table 5: Values for c = dv{H*)/dH*\H*-o, where v{H*) is the volume fractfon 
of one of the two labyrinths for the corresponding family of surfaces of constant 
mean curvature. The measured values follow from our sphne interpolation of 
the numerical data of Ref. [2]. In the second row we give —Ao'^/2'kx, a new 
estimate for c derived in the text. 

H = Cq minimize the functional / dA{H — cq)'^ . In any case, since H is constant 
over the surface, surface area and volume now vary in the same way, and we 
have dA/dV = 2H. 

All simple TPMS which are of interest for physical reasons are members of 
a family of CMC-surfaces [55]. Each of these families consists of two branches, 
corresponding to positive and negative mean curvatures, which are separated 
by the minimal surface member. Like the minimal surface member, each of the 
the triply periodic CMC-surfaces of the family partitions space into two inter- 
twined, yet separate labyrinths, with volume fractions v and I — v. Each family 
of CMC-surfaces, and therefore the volume fractions of the two labyrinths and 
the surface area A, is parametrized by H. If the TPMS-member is balanced, 
then the operation a maps one branch of the family onto the other, in partic- 
ular v{H) = 1 - v(-H), v{H = 0) =vo = 0.5 and A{H) = A{-H). For an 
amphiphilic monolayer of thickness (5 in a ternary system, v can be identified 
with the hydrocarbon volume fraction (that is oil and amphiphilic tails), 1 — v 
with the water volume fraction (where we neglect the contributions of the am- 
phiphilic heads), and 2A5 with the amphiphilc volume fraction. In the seminal 
work by Anderson, Davis, Nitsche and Sciven from 1990, the CMC-families were 
numerically constructed for P, D, I-WP, F-RD and C(P) using finite element 
methods [2]. In 1997, Grosse-Brauckmann numerically constructed the G-family 
in a similar way (using the software package Surface Evolver) [44]. For these 
families, the following data has been tabulated: volume fraction ?; as a function 
of scaled mean curvature H* and scaled surface area A* as a function of v. 
Rearrangement and interpolation with cubic splines provides smooth functions 
H*{v) and A*{v). As v is varied away from the value vq for the TPMS (= 0.5 
for balanced structures), one has 

H*{v) = -^^^^^ + 0{{v-v^f) (23) 
c 

A*{v) = Ao- + 0{{v- vof) 

c 

where the numerical value of c can be extracted from the cubic splines (compare 
Tab. 5). Note that the two relationships are not independent due to the relation 
dA/dV = 2H. They are related through c, whose values are well approximated 
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by the analogous values for the parallel surface companions to the TPMS, which 
can be derived as follows [98] : if 6 denotes the perpendicular distance from the 
minimal to its parallel surface, to lowest order in 6 the volume fraction v and 
the mean curvature H* , averaged over the surface in the unit cell, are given by 
V = vo + A*S and H* = 2TrxS/A*, respectively. Thus H* = 2Trx{v - vo)/A*^ 
and c = —A*'^/2'KX for the parallel surface case. The corresponding numbers 
arc given in Tab. 5; except for C(P), the overall agreement with the numerical 
data for c for the CMC-surfaces is remarkably good. 

Eq. (23) is a useful approximation for CMC-surfaces close to the TPMS, 
where they behave similar to parallel surfaces. However, as mean ciirvature 
grows, the numerical data starts to deviate from these approximations, changes 
in V and A* become slower, and a turning point is reached, where v reaches 
an extremal value and starts to decrease again as a function of H. Beyond the 
turning point, the surfaces correspond to nearly spherical regions connected by 
small necks which resemble pieces of unduloids. Finally these necks disappear 
and each branch terminates in an assembly of sphere, which might be close- 
packed or self-intersecting. Fig. 9 shows the geometrical data for the families G, 
D, P, C(P), I-WP, F-RD as tabulated in the literature [2, 44]. We do not show 
the parts of the data beyond the turning point, as these correspond to surfaces 
whose structure is unphysical. 

We now consider a ternary mixture of water, oil and amphiphile, where am- 
phiphiles self-assemble into monolayers with spontaneous curvature cq . In such 
a system, all relevant phases (micellar, hexagonal, lamellar, cubic bicontinuous) 
can be modeled as CMC-surfaces (spheres, cylinders, planes, triply periodic 
CMC-surfaces). This approach has first been used by Safran and coworkers 
[94, 95, 116], who in particular discussed the surfaces from the D-family. In our 
work [98], we extended this analysis to all families of interest, including the G- 
family, for which the relevant data has become available only recently and which 
features very prominently in experimental systems. A ternary mixture has two 
independent degrees of freedom for concentration, which we choose to be the 
hydrocarbon volume fraction v and the ratio w of amphiphile to hydrocarbon 
volume fraction. For the formulae given later it is in fact useful to modify the 
definition of w and to scale it with the dimensionless spontaneous curvature cq: 
w = pa/vcqS. The analysis of Eq. (1) for CMC-surfaces is simplified consider- 
ably by the fact that as mean curvature H is constant, there is no difference 
between local and global curvature properties and the integral over the surface 
becomes trivial. We use a factor 2kCq to rewrite the free-energy density in the 
dimensionless form 



For phases with lamellar, cylindrical and spherical aggregates, this free-energy 
density can easily be expressed as a function of the concentration degrees of 
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freedom: 
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Note that only fs depends on r, since the other two structures have no Gaussian 
curvature. Since the aggregates are disconnected, the dependence on v is trivial. 
The phase boundaries S — C, S — L, C — L and the emulsification failure are 
obtained from Eqs. (25), (26) and (27) to be w = 24/(7 - 16r), w = 6/(1 - r), 
w = 8 and w = 3/(r — 1), respectively. For example, for r = and increasing 
w, we find the phase sequence L - C - S - emulsification failure, which is typical 
for amphiphilic systems (for simplicity, we identify phase transitions with cross- 
ing points of the free-energy-density curves, and use the Maxwell construction 
only for the emulsification failure). In Fig. 4 one sees that experimentally the 
emulsification failure indeed occurs at constant w (straight line through water 
apex). For increasing r, the spherical phase becomes more favorable and finally 
suppresses the cylindrical phase. 

The structures based on CMC-surfaces discussed here are very different from 
the lipid bilayer structures discussed in the preceding section: now only one 
monolayer is present, and one of the two labyrinths is filled with hydrocarbon. 
In the following, they are called single structures. It should be noted that 
for non-balanced structures, it makes a difference which of the two labyrinth 
is filled with hydrocarbon; for example, the I-WP-family generates two single 
structures, which we call I and WP. In ternary amphiphilic systems there also 
exists an analogue to the bilayer structures discussed before, which we call 
double structures and mark with an index /. Double structures have the same 
geometries like cubic bicontinuous phases in diblock copolymer systems. They 
can be considered to be TPMS-based bilayer structures where the inner part 
of the bilayer has been swollen with suitable solvent. Since each of the two 
monolayers has the same spontaneous curvature cq, a double structure can be 
modeled as the combination of the two surfaces of a CMC-family with H = Cq 
and H = —cq. Since for not too large H these two CMC-surfaces essentially 
correspond to the shrinkage of one of the two separate labyrinths defined by 
the minimal surface member, the two interfaces of a double structure do not 
intersect. Note that in principle there exists another class of structures, that 
is double structures where oil and water (and therefore amphiphile orientation 
and mean curvature) have been reversed. However, since finite spontaneous 
curvature selects only phases with curvature towards one specific side, they are 
suppressed for energetic reasons by the lamellar phase. 



In our work [98], we considered 8 different single structures, which exist for 
the volume intervals [0.056,0.944] for G, [0.131,0.869] for D, [0.249,0.751] for 
P, [0.481,0.519] for C(P), [0.357,0.857] for I, [0.143,0.643] for WP, [0.439,0.625] 
for F and [0.375, 0.561] for RD. We also considered 6 double structures, which 
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exist for the volume intervals [0.112,1.0] for G/, [0.262,1.0] for D/, [0.498,1.0] 
for Pi, [0.962,1.0] for C(P)/, [0.624,1.0] for I-WP/ and [0.818,1.0] for F-RD/. 
Note that the gyroid structures cover the largest intervals in v for their respec- 
tive class: there is no other structure which can incorporate so extreme volume 
fractions like the gyroid. Together with the 3 non-cubic phases treated above, 
we considered 17 different phases. Since the cubic phases consist of one con- 
nected aggregate, the scaling of the free energy density / with hydrocarbon 
volume fraction v will be non-trivial. For a given value of v, the mean curvature 
H{v,a) = H*{v)/a and the surface area A{v,a) = A*{v)a'^ within a unit cell 
are determined by the curves plotted in Fig. 9. The amphiphile concentration 
pA = A{v,a)/a^ = A*{v)/a fixes a, so that a = A*{v)/pA = A* {v) / {wvcq) . 
Prom Eq. (24) we can then derive 

fsciw, V, r) = wv [K{v) wv-l] + r (28) 

where we have used the definition for the curvature index A and the topology 
index F from Eq. (7). Note that now the indices are f -dependent. For the 
single structures, / follows by combining using the data shown in Fig. 9 with 
the definitions in Eq. (7). For the double structures, the procedure is somehow 
more complicated. However, for balanced double structures, it becomes simple 
again, one then can use 

Ki{v) = AK2)/2, r/(^;) = 2r(^;/2), fscAw, v, r) = 2fBc{w, ^, r) (29) 

in Eq. (28). 

For r = (vanishing saddle splay modulus R), the analysis of cubic bicon- 
tinuous phase behavior becomes rather simple: in contrast to the case of the 
parallel surface model discussed in the last section, now the bending energy given 
in Eq. (28) can be completely relaxed, namely by satisfying w{v) = l/wA(w). 
This leads to lines of vanishing frustration in the Gibbs triangle. For the free- 
energy density of the micellar and hexagonal phases (Eq. (27) and Eq. (26), 
respectively), the lines of vanishing frustration follow as w = 3 and w = A, 
respectively. All lines of vanishing frustration are plotted in Fig. 10a. Obvi- 
ously for this case phase behavior is very complex and degenerated, since every 
structure considered has some region of stability around its line of vanishing 
frustration. This type of degeneracy caused by the bending energy has been 
discussed before [6] and leads to the conclusion that additional physical effects 
have to be operative, as it is not observed in experimental systems. 

Using Eq. (23), a simple approximation can be derived for the lines of van- 
ishing frustration of the cubic phases: 

w = -cAo/va^v - vo), wj = -4c^o/(^' - 1) • (30) 

Thus the minimal surface case corresponds to the stable solutions for w ^ 1 
at V = vo and v = 1, respectively. For the balanced single structures and 
the double structures the hierarchy of the different phases within the band-like 
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(a) 



(b) 



Figure 10: (a) Lines of vanishing frustration for r = 0. With increasing am- 
phiphile concentration, the structures S - C- double cubic - single cubic are 
stable, (b) Phase behavior for r — 1/15. Now all cubic phases have been 
suppressed except the two gyroid phases. In (a) and (b), cq = 1/6, which 
corresponds to H2O /Cii/Ci2E5 at T = 20°C. 

region occupied by a certain structural type is thus determined by the values 
of cAq. Using the approximation c ~ — A*^/27rx derived above, we find cA* « 
r^, where A* and T correspond to the minimal surface members. Therefore 
the phase sequence is approximately determined by the topology index of the 
minimal-surface member of each family. In particular, for a given structural 
type we expect to find the sequence G - D - P as a function of either v or w. 

For r > (negative saddle splay modulus R) , spheres become more favorable, 
cubic phases less favorable and cylinders and lamellae experience no change in 
free energy. Therefore the cubic phases will finally disappear, but our numer- 
ical analysis shows that cubic phases can persist up to r = 0.2, in contrast 
to earlier work, which predicted r = 0.1 [116]. The reason for this becomes 
clear in Fig. 10b where we show the full phase diagram for r = 1/15 = 0.07: 
the only cubic phases stable here are the two gyroid structures, which have 
not been considered before. There are two main reasons for their outstanding 
performance: Eq. (28) shows that their large values for the topology index F 
reduces the energetic penalty caused by r, and since they can accommodate ex- 
treme volume fractions, they can compete with other structures at all relevant 
concentrations. By comparing the theoretical phase diagram from Fig. 10 with 
the experimental one from Fig. 5, we conclude that the experimental system 
should correspond to a rather large value of r. Earlier work moreover suggests 
that incorporation of thermal fluctuations would favor micellar phases near the 
water-amphiphilic side and lamellar phases near the water apex [101]. Such a 
modification is expected to considerably improve the agreement between the 
two phase diagrams. 
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Figure 11: Freeze- fracture microscopy picture of a balanced microemulsion 
phase. From Rcf. [52]. 

6 Random surfaces 

6.1 Microemulsion and Sponge Phases 

When K/fc^r, the bending rigidity in thermal units, becomes sufficiently low 
(that is temperature T sufficiently high), bicontinuous structures can now longer 
maintain their long-range crystalline order and melt into a disordered phase, 
which is characterized by an exponential decay of correlations in the interfa- 
cial positions. Such phases have been observed experimentally for a long time 
in many binary and ternary amphiphilic systems. Microemulsions are macro- 
scopically homogeneous and optically isotropic mixtures of oil, water and am- 
phiphiles. On a mesoscopic scale, they consist of two multiply connected and 
intertwined networks of oil- and water-channels, which are separated by an 
amphiphilic monolayers. Free-fracture microscopy, where the sample is quickly 
frozen, cut, and then studied with an electron microscope, reveals the intriguing 
structure of this phase [52], see Fig. 11. A similar phase, the sponge phase, ap- 
pears in binary systems of water and amphiphile, where now the two labyrinths 
are occupied by water, which are separated by an amphiphilic bilayer. The pic- 
tures obtained by free- fracture microscopy [107] are even more suggestive in this 
case, because the sample has a preference to break along the bilayer mid-surface, 
so that the three-dimensional structure of the membrane becomes visible. An 
example is shown in Fig. 12, which clearly shows the saddle- like geometry of the 
amphiphile film. Therefore, the intuitive picture of microemulsion and sponge 
phases as fluid versions of bicontinuous cubic phases is strongly supported by 
these experiments. 

These phases have been investigated experimentally in considerable detail 
over many years. In particular, their phase behavior and scattering intensi- 
ties have been studied carefully. A theoretical understanding of the statistical 
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Figure 12: Freeze-fracture microscopy picture of a sponge phase. From 



mechanics of membranes, however, is only beginning to emerge in recent years 
[75, 40, 78, 18, 35]. This is no surprise, since the statistical mechanics of a sur- 
face, which can not only change its shape, but also its topology in all possible 
ways, is extremely complicated. In principle, a partition function of the form 



has to be calculated, where I?R(r) denotes an integration over all possible 
shapes with parametrization R(t) of the surface at fixed topology, where r is a 
two-dimensional coordinate system on the surface. However, this integral cannot 
be just over all possible parametrization R(t) of a surface of fixed topology, but 
has to be restricted to those parametrizations, which lead to physically different 
shapes in the embedding space; this is indicated by the prime. Finally, the 
contributions off all different topologies have to be summed over. It is clear 
that this problem is sufficiently complex that no exact solution will be found 
anytime soon. Therefore, approximations have to be made in order to get some 
insight into the behavior of these phases. 

6.2 Gaussian random fields 

A very useful approach is to describe the interfaces as isosurfaces of Gaussian 
random fields (GRF). This corresponds to a Ginzburg-Landau model as dis- 
cussed in Sec. 2.2, in which the free-energy functional is taken to be quadratic 
in the scalar field <I>(r). This Gibbs distribution is homogeneous and isotropic, 
and all the functional integrals are Gaussian and can be performed exactly. In 
fact, there is a considerable mathematical literature on the isosurfaces of GRF, 
see e.g. the book by Adler [1]. In order to comply with the conventions of this 
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field, wc now make two changes to our notation. In tlie following, area con- 
tent is denoted by S (rather than by A) and amphiphilic (hydrocarbon) volume 
fraction by ^ (rather than by v). 

The amphiphilic monolayers in microemulsions have been modelled as level 
surfaces of GRF by Berk [4, 5], Teubner [112], Pieruschka and Marcelja [81], and 
Pieruschka and Safran [82, 83]. Here the starting point is a Gaussian free-energy 
functional of the general form 



Ho[</.] 



i I ci'^gz/(q)-i$(q)<l>(-q) . 



(32) 



The average geometry of the $(q) = a level surfaces can be calculated for arbi- 
trary spectral density i'(q). For the surface density, S/V, the mean curvature H, 
the Gaussian curvature K, and the mean curvature squared H'^, the following 
averages are obtained [112]: 
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The value a of the level cut can be used to describe the preferred curvature of 
the membrane, as well as the volume fractions of oil and water. Since {H) is 
a linear function of a, compare Eq. (35), this parameter is proportional to the 
spontaneous curvature cq. In particular, for a = the mean curvature of the 
surface vanishes; this applies to a balanced system, where cq = 0. 

For balanced systems, it follows from Eqs. (33,34) that the topology index 
r = V^/tt = 0.9003 (compare Eq. (7)). This value is only slightly larger than 
the ones for cubic bicontinuous phases (compare Tab. 1), since the balanced 
random sponge features only few disconnected parts (note that F doubles when 
the structure is duplicated). For small curvatures (that is small a), one can use 
Eqs. (33,34,35) to derive a relationship between topology index F and curvature 
index A which is independent of spectral density !/(q) and a: 



(38) 



When the curvature index A increases since the sponge's interfaces gain curva- 
ture, the topology index increases, too, since disconnected parts proliferate. 
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The GRF-approach is most predictive when the Gaussian model of random 
interfaces is related to the statistical mechanics of membranes by a variational 
approximation [82, 83]. In this case, the spectral density u{q} in the functional 
(32) is determined by the requirement that the (f){r) = level surfaces mimic 
the behavior of interfaces controlled by the curvature Hamiltonian (1) as close 
as possible. The usual variational approach employs the Feynman-Bogoljubov 
inequality, 

F <Fo + {n-Ho)o (39) 

where 7i and F are the Hamiltonian and the free energy of the system of in- 
terest, respectively, and Ho and Fq the same quantities of the reference system. 
In this way. an upper bound for the true free energy is obtained. This is more 
complicated in the case of random surfaces, because the GRF-Hamiltonian is 
defined everywhere in space, while the curvature Hamiltonian is only defined 
on the level surface. Therefore, the curvature energy does not restrict fluctu- 
ations of the field $(r) away from the level surface. In order to suppress such 
fluctuations, one usually makes the mean-spherical approximation, that is the 
constraint ($(r)^) = 1 is introduced. 

With this variational approach, Pieruschka and Safran [82] have been able 
to derive the following form for the spectral density 

and to relate the parameters a, b and c to the curvature elastic moduli k and 
K and the surface density S/V. The first interesting result is that the spectral 
density is found to be independent of R. Exact expressions for the parameters 
can be found in Ref. [22]. To leading order in ksT/n, the parameters simplify 
to 

15^2 S 
a = , (41) 
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The spectral density (40) is equivalent to the scattering intensity in bulk con- 
trast. Its Fourier transform yields the correlation function 

($(r)$(r')) = J d^q e*'i<'-'''^u{ci) = j exp[-r/^] sin(fcr) (44) 

where the second equality holds for |6| < Thus, the correlation function 

is characterized by two length scales, the correlation length ^ and the typical 
domain size d = 27r/fc of the oil- of water channels, which are obtained from 
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Eqs. (40) and (44) to be [113] 



= ^pV~c + h, (45) 



k 



= y^Vc-b. (46) 



With the results (41), (42) and (43), the asymptotic behavior for small ksT/n 
of the dimensionless product k£^ is found to be 

The free energy of the sponge phase can also be calculated from the GRF ap- 
proach. To leading order in an expansion in ksT/K, the free-energy density 
/ = F/V is found to be [83] 

/^^i— i(f)'-¥K4)- 

This implies that for small membrane volume fractions = 5S/V (where again 
5 is the thickness of the amphiphilic interface), the entropic term dominates over 
the energy term, and that the sponge phase becomes unstable in this regime. 

Wc want to end this section with a short discussion of the reliability of the 
predictions of the GRF model as a variational approximation for membrane 
ensembles. The weak point of this approach is that it is not clear whether the 
calculated entropy is actually equivalent to the physical conformational entropy 
of the membranes [74]. The main problem is that the curvature energy only 
controls the shape of the $(r) =0 level surface, while the values of the scalar 
field $ at all other points in space are not affected by it. The fluctuations of 
$ in these oil- and water-regions are mainly determined by the mean-spherical 
constraint ($^(r)) = 1. Obviously, an appreciable contribution to the total 
entropy arises from the fluctuations of <I> these 'bulk' regions. This would not 
affect the predictions of the model as long as the 'bulk' contributions were 
independent of the interface positions. Unfortunately, there is no argument so 
far that this is indeed the case. 



6.3 Phase behavior of random surfaces 

From here on we consider only the case of vanishing spontaneous curvature, 
Co = 0, that is balanced microemulsions and sponge phases. Then the phase 
behavior is controlled by the bending rigidity k, the saddle-splay modulus k, 
the density of membrane area per volume, S/V , and a microscopic cutoff (5, 
which can be identified with the thickness of the amphiphilic interface. It is 
very important to realize that for vanishing spontaneous curvature, phase tran- 
sition as a function of the amphiphile volume fraction cannot be understood on 
the basis of the curvature energy alone — i.e. without considering the effect 
of thermal fluctuations [85]. The reason is that the curvature Hamiltonian is 
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confornially invariant in three spatial dimensions, which imphes in particular 
that it is invariant under a simultaneous rescaling of all length scales. Since 
the curvature energy is scale invariant, the energy density scales as the third 
power of an inverse length, i.e. as (S/V)^. Therefore, the curvature energy 
of any given structure — spherical, cylindrical, lamellar, cubic, or random — 
scales in exactly the same way with decreasing amphiphile volume fraction, and 
their relative order is maintained. Therefore, thermal fluctuations are crucial 
for these phase transitions. 

It has been suggested by several authors [93, 9, 73, 28] that the free energy of 
the sponge phase can be obtained by integrating out the membrane fluctuations 
on scales less than the typical domain size. This integration over small-scale 
fluctuations leads to renormalized, scale-dependent curvature moduli Kn{l) and 
kr{1) as given by Eqs. (4) and (5), respectively, with l/S = where ^' is 

the membrane volume fraction. This implies that the curvature Hamiltonian 
(2) has to be replaced by 



F = J dA|iK+,fl(vE'-l)(fci+fc2)' + ^/^-,ii(*-')(/Cl-fc2)'| 



(49) 



with 



K+,fl(M/-i) = «+ + M^ln*, (50) 

OTT 

/c_,fl(*-i) = K_ + ^ln*. (51) 

The stability arguments used in Sec. 2.1 imply that k+^r and k^,r have to be 
positive for the free energy (49) to be stable against collapse of the structure 
to molecular scales. Therefore, there are instabilities at K+^i{(\J'~^) = and 
K-^r{'^~^) = 0. The latter instability can be identified with the emulsification 
failure of the sponge phase, so that the phase boundary is predicted to occur at 
[73] 

, -r 67r K 

This resiilt can be understood intuitively as follows. For sufficiently large mem- 
brane volume fraction, both and K-^r are positive. Therefore, the system 
tries to minimize (fci + ^2)^ and (ki — ^2)^. This can be achieved by decreasing 
both ki and ^2, i.e. by swelling a given structure as much as possible the 
lamellar or sponge phase is stable at this value of On the other hand, as soon 
as K;+,ij or K-^r become negative at some small value of ^, the free energy can 
be reduced by collapsing the structure. With decreasing length, however, k+,_r 
and K-^R increase and finally become positive. Therefore, the collapse stops at 
a length scale which is exactly the length scale determined by Eq. (52). 

It is worth mentioning that a similar result follows from the calculation of 
the shapes and free energies of passages in lamellar phases [31]. Passages are 
catenoid-like connections between adjacent lamellae which proliferate close to 
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phase transitions to disordered bicontinuous phases. From a detailed calcula- 
tion, which takes into account the Gaussian membrane fluctuations, the density, 
p, of passages per unit base area of the stack can be obtained to be 
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Here, A(j(K/kBT) is a function which has to be calculated numerically; it is pre- 
dicted to increase monotonically with decreasing n/ksT in Ref. [31]. Therefore, 
the density of passages increases with increasing lamellar spacing d/5 = *~^. 
When the average distance between passages shrinks to the average membrane 
separation, a transition to the sponge phase can be expected to occur. This 
happens at 



where ** is a fimction of n/ksT. Equations (52) and (54) agree perfectly. 

The consideration of the renormalization of k and k implies that the free- 
energy density, /, of the sponge phase should behave as [89, 87, 90] 



As explained above, the overall scaling with derives from the conformal 
invariance of the curvature energy. A represents the bending energy without 
thermal fluctuations, and is a linear function of both k and R. B represents 
the logarithmic corrections from the renormalization and is a linear function of 
temperature. Both A and B depend on the detailed geometrical structure of 
the sponge phase, and thus cannot be obtained from simple scaling arguments. 
It is important to note that the functional dependence of the free energy (55), 
which is based on the renormalization of the curvature elastic moduli due to 
small-scale membrane fluctuations, does not agree with the free energy (48) of 
the Gaussian random field model, which includes the topological entropy of a 
disordered bicontinuous phase. 

6.4 Monte Ceirlo simulations of triangulated surfaces 

In order to go beyond the approximations discussed in Sees. 6.2 and 6.3, dis- 
cretized surface models can be investigated by Monte Carlo simulations [35, 34, 
36, 37]. Surface triangulations provide the best way of discretizing a surface as 
uniformly as possible. The model consists of vertices, which are connected by 
bonds in such a way that the bonds form a triangular network. Several Hamilto- 
nians have been suggested for triangulated surfaces, such that their shapes and 
fluctuations are governed in the continuum limit by the curvature energy from 
Eq. (1) [35]. Two examples are Gaussian-spring models, in which neighboring 
vertices — the vertices connected by bonds — interact with harmonic spring 
potentials, and tether-and-bead models, in which the interaction V{r) between 
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neighboring vertices is defined the potential 

{oo < r < (To 

for CTo < r < io (56) 
00 ^0 <r 

where ctq is the hard-core diameter of the beads, and (.q is the tether length. For 
tether lengths to < \/3cro, the network is self- avoiding, because a bead of some 
distant part of the membrane does not fit through the hole between the three 
beads of a triangle, even for maximally stretched tethers. 

A Monte Carlo simulation of sponge phases requires three types of Monte 
Carlo movers [36] . The first step is to displace selected beads by a random vector 
chosen uniformly in the cube [— s, s]"^. Here, s is determined by the criterion 
that roughly 50% of trial moves is accepted according to the Boltzmann weight. 
The second step is to vary the connectivity of the bond network, in order to 
allow for diffusion and fluidity of the membrane. This dynamic triangulation is 
performed for any two adjacent triangles. The bond which forms the common 
edge of the two triangles is cut, and a new bond is inserted, which connects 
the two previously unconnected beads. The bond can only be cut if every 
bead retains bonds to at least three other beads. This procedure guarantees 
that the network remains two-dimensionally connected, no holes open up in the 
network, and the topology of the network does not change. Finally, the third 
step is to change the topology of the surface. This is done by removing two 
nearby triangles from the surface, and by connecting the corresponding vertices 
by a prism of six new triangles. Of course, the inverse step is also possible, 
with a passage six triangles being removed, and two new triangles inserted to 
close the surfaces. The acceptance of both the bond-flip and topology-change 
moves is determined by the Boltzmann weight and is therefore controlled by the 
curvature energy. We want to mention parenthetically that some care has to be 
taken to find a good discretization of the bending energy [33] . 

A typical configuration of a triangulated surface in a cubic box with param- 
eters K, K and \I' chosen in the stability region of the microemulsion or sponge 
phase is shown in Fig. 13. This configuration nicely demonstrates the bicontin- 
uous structure of balanced microemulsions and sponge phases. The saddle-like 
geometry of the membrane can also be easily seen. Finally, the figure shows that 
locally the structure of the sponge phase strongly resembles the cubic phases 
discussed above. Therefore, a sponge phase should indeed be considered as the 
molten state of the crystalline cubic phase. 

A more quantitative comparison with the theoretical approaches of Sees. 6.2 
and 6.3 above can be made by determining the phase diagram of the randomly- 
triangulated surface model, and by calculating the osmotic pressure, p, in the 
simulations as a function of the membrane volume fraction From Eq. (55), 
we obtain 

= [{2A + B) + 2B\n^^^ (57) 
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Figure 13: A typical membrane configuration in a sponge phase for bending 
rigidity K/ksT ~ 1.6. The two sides of the membrane are shaded differently in 
order to emphasize the bicontinuous structure of this phase. From Ref. [36]. 

i.e. the same functional dependence as the free-energy density itself. This de- 
pendence of the osmotic pressure is indeed nicely confirmed by the simulation 
data [36]. The simulations therefore provide strong evidence for the renormal- 
ization of the elastic moduli of the curvature model and for the dependence (55) 
of the free energy on the membrane volume fraction. 

The phase diagram for fixed bending rigidity k is shown as a function of 
K and 5" in Fig. 14. The simulation data are compared with the prediction 
(52) for the phase boundary. Since the slopes of the phase boundaries in this 
logarithmic plot agree very well, not only the exponential dependence of the 
membrane volume fraction at the transition on the saddle splay modulus is 
confirmed, but also the value of the imiversal prefactor in Eq. (52) is strongly 
supported. 

6.5 Comparison with experiments 

Experimentally, phase diagrams and scattering intensities have been studied 
systematically for many different surfactant molecules. Qualitatively, the agree- 
ment with the theoretical approaches is very reasonable. For example, the 
scattering curve in bulk contrast shows a peak at non-zero wave vector in the 
microemulsion phase, which moves out and decreases in height with increasing 
surfactant concentration. A quantitative comparison, however, is much more 
difficult. We want to discuss here three different classes of experiments, where 
a such a quantitative comparison has been made. 

The first type of experiments are scattering studies in bulk and film contrast. 
In bulk contrast, the scattering intensity Swwi<i) is proportional to the spectral 
density i/(q) of the Gaussian random field model, compare Eq. (40), for wave 
vectors q which are not much larger than the characteristic wave vector k of 
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Figure 14: The phase diagram as a function of membrane volume fraction <J) and 
saddle-splay modulus R, for bending rigidity k/UbT ~ 1.6. Note the logarithmic 
scale of the abscissa. The dashed line shows the theoretical prediction (52). 
Prom Ref. [36]. 

the domain structure, compare Eq. (46). The functional dependence of Eq. (40) 
describes the scattering data in this regime very well, as was first noted by 
Teubner and Strey [113], who derived this result on the basis of a Ginzburg- 
Landau model, very similar to that introduced in Sec. 2.2. We want to mention 
parenthetically, that for wave vectors k ^ q <^ 1/5, the intensity is dominated 
by the scattering from sharp, planar interfaces; in this limit, the famous Porod 
law [84, 112] predicts I{q) - {S/V)q-'^. 

In the limit of wave vector q ^ 0, the scattering intensity in film contrast is 
given by [87] 

' (58) 

where p is the osmotic pressure of Eq. (57). For the free energy (55), this implies 

[*J(g ^ 0)]-^ ~ const -h In * (59) 

Such a behavior has indeed be observed experimentally in Ref. [87]. However, 
this result has been questioned by Daicic et al. [14, 13]. This has lead to a 
intensive debate, with arguments against [15, 16] and in favor [91, 86] of the 
existence of a logarithmic renormalization of clastic moduli in the sponge phase. 

In a second type of experiment, information about the average geometry of 
the surfactant film can be extracted from the scattering intensity in the regime 
k < q 1/5. This information is contained in the corrections to the asymptotic 
1/g^ law for smaller values of the wave vector [111]. Experimentally, the average 



38 



Gaussian curvature is found to be [11] 

{V/Sf{K) = l.2f,±Q.lQ (60) 

which is in excellent agreement with the Gaussian random field result (V/ S)"^ {K) = 
1/r^ = — 7r^/8 = 1.23 and simulations of the $^-Ginzburg-Landau model 
[32, 30]. 

The third type of experiment concerns the phase behavior of mixtures of 

water and non-ionic surfactant as a function of temperature and surfactant 
concentration. As mentioned in Sec. 4, the saddle-splay modulus k is a linear 
function of temperature in this case. Therefore, the concentrations at the phase 
boimdaries of the lamellar and the sponge phase are expected from Eq. (52) 
to depend exponentially on temperature. A logarithmic plot of the phase di- 
agram of C12-B5 in water [109], compare Fig. 4, is indeed consistent with this 
expectation. 

The most detailed information about both the scattering intensities and 
phase behavior in these systems has been obtained very recently in ternary 
amphiphilic systems of water, oil and non-ionic surfactant CiEj, to which small 
traces of an amphiphihc block copolymer has been added [53, 21, 22, 29, 38]. 
The results obtained in this system provide further, strong evidence for the 
existence of a logarithmic renormalization of k and R in microemulsions and 
sponge phases. 

7 Summary and outlook 

In this contribution we discussed the geometrical properties of surfaces which 
can be used as structural models for cubic bicontinuous phases in amphiphilic 
systems: triply periodic minimal surfaces, their parallel surfaces and constant 
mean curvature companions, and bicontinuous random surfaces. For each class 
of surfaces, we showed how the geometrical properties translate into physical 
properties of bicontinuous phases in amphiphilic systems in the framework of 
an interfacial description. The surprising success of this approach relates to the 
fact that in amphiphilic systems, the solvent has little physical properties by 
itself and the free energy is essentially determined by the interfaces. Although 
there are several physical effects which have been neglected in our treatment of 
amphiphilic systems, including van der Waals and electrostatic interactions, it 
can be concluded that the most essential aspects of phase behavior are now well 
understood in terms of the properties of the underlying geometries. 

During the last years, the interest in bicontinuous phases has increased due 
to some promising applications in the nano- and biosciences. For example, am- 
phiphilic self-assembly has been used to synthesize mesoporous systems [58], 
that is porous material with amorphous walls (usually silica-based) and pore 
sizes on the nanometer scale. As a matter of fact, the combination of am- 
phiphilic self-assembly and crystallization also seems to be a basic aspect of 
biomineralization [70], and there are many algae whose mineral skeletons look 
similar to cubic bicontinuous phases [51]. There is also a large effort underway to 
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synthesize bicontinuous structures from graphitic material (Schwarzites), which 
offers the advantage of atomically smooth walls [69]. Porous material on the 
nanometer scale is not only interesting for its structural properties (for example 
use as sieves or catalysts) , but also for its electronic properties (for example use 
as photonic bandgap material). 

In the biosciences, two recent developments involved cubic bicontinuous 
phases. It was found that the gyroid phase in the monoolcin-watcr system 
provides a functional environment for the crystallization of integral membrane 
proteins like bacteriorhodopsin [77, 66]. Membrane proteins are notoriously 
difficult to crystallize in a three-dimensional array, which however is necessary 
for structure determination by X-ray scattering. And an extensive analysis of 
transmission electron micrographs of many biological specimen has shown that 
cubic bicontinuous structures are locally formed in several regions of the cells 
[61, 19], including endoplasmatic reticulum, Golgi apparatus and mitochondria. 
These are extended lipid bilayer systems in the cell which carry a protein ma- 
chinery responsible for processing and packaging material inside the cell. In fact 
it is easy to imagine that biological cells regulate their lipid composition and 
therefore their spontaneous curvature as to achieve geometries favorable for the 
task at hand. In general, cubic bicontinuous phases seem to be ideal space par- 
titioners if a large amount of active surface area is needed while simultaneously 
providing good access for incoming and outgoing material through the space 
away from the interface. 

What problems are left regarding the subjects of this article ? From the 
viewpoint of mathematics, there are two interesting aspects. On the one hand, 
the physical motivation leads to certain variational problems for which little 
is known about the corresponding solutions. One example is the energy func- 
tional J dA{H — co)^ with a volume constraint; in principle this problem was 
subject of Sec. 5, but in lack of adequate mathematical knowledge (this prob- 
lem has been treated for vesicles [104], but not for bicontinuous structures), 
we used CMC-surfaces as structural models. On the other hand, some classes 
of surfaces are well investigated, but still representations are missing which 
could be used in physical modelling. Therefore additional progress in regard 
to Weierstrass representations for TPMS, Weierstrass-like representations for 
their CMC-companions, and other (numerical) representations would be most 
welcome. With a larger repertoire of representations, it should also become 
easier to treat non-local effects like van der Waals and electrostatic interac- 
tions; here integral geometry might be helpful in reducing the dimensions of the 
corresponding integrals, and sophisticated summation techniques are needed to 
account for the infinite extension of the bulk phases. 

From the viewpoint of material sciences, it seems worth to study the role 
of boundaries and inclusions. The term boundaries includes grain boundaries 
between different, yet adjacent phases, for example between a hexagonal and 
a gyroid phase; until now, this problem has been investigated only for simple 
geometries. It also includes the boundary to an external surface (which might 
even be curved or chemically structured) and free boundaries. For example, it 
has been recently observed experimentally that the bicontinuous cubic phase 



40 



in the binary system water - C12EQ is the only known material system which 
shows crystal faces with high Miller indices {devil's staircase) [80]. The term 
inclusions includes several object of biological and technological relevance, in 
particular integral membrane proteins, polymers anchored in the membrane and 
(functionalized) nanoparticles. In all of these cases, the framework described 
here should provide a good starting point for future research, which has to bridge 
the gap between the mathematical knowledge on the geometrical properties and 
the physics knowledge of the material basis of these fascinating structures. 

8 WWW resources 

Here we list some internet addresses which relate to the subject of this contri- 
bution: 

• http : //www-sf b288 . math . tu-berlin . de/~konrad/ : 

Homepage of Konrad Polthier. one of the pioneers for numerical represen- 
tations and visualizations of minimal surfaces. 

• http : //people . math . iini-bonn . de/kgb/: 

Homepage of Karsten Grosse-Brauckmann, a mathematician working on 
CMC-surfaces, in particular on the G-family. 

• http : //www . susqu . edu/f acstaf f /b/brakke/evolver/: 

The Surface Evolver is a public domain software package written by Ken 
Brakke for the propagation of triangulated surfaces. In particular. Ken 
Brakke shows many examples of TPMS obtained with the Surface Evolver. 

• http : //www.msri . org/publications/ sgp/SGP/indexc .html: 

The Scientific Graphics Project oH'ers stunning visualizations of minimal 
and CMC-surfaces, obtained mostly in the framework of Weierstrass rep- 
resentations. 

• http://www.gajig.umass.edu/: 

The Center for Geometry, Analysis, Numerics and Graphics (GANG), 
where a lot of work on CMC-surfaces takes place. 

• http : / / www . mpikg-golm .mpg . de/th/people/ schwarz/ : 

More information on the work described here, including color pictures of 
all structures investigated and (improved) nodal approximations. 
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